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§ Abstract 

O The question of decoupling and freeze-out is reinvestigated and analysed in terms 

^ of transparent semi-classical decoupling formulae, which provide a smooth decou- 

pling in time both, for single and two particle inclusive spectra. They generalise 
frequently employed instantaneous frcczc-out procedures and provide simple rela- 
tions between the damping width and the duration of the decoupling process. The 
^ implications on physical phenomena arising from the expansion and decay dynamics 

^ of the highly compressed hadronic matter generated in high energy nuclear collisions 

are discussed. 
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O 1 Introduction 

> 

^ Dynamically expanding systems such as those created in high energy nuclear 

collisions or the early universe pass various stages, where different degrees of 
freedom are relevant and thus different dynamical concepts are appropriate. In 
the context of nuclear collisions, cf. [1], one envisages a phase, where at high 
densities partonic degrees of freedom prevail in form of a quark gluon plasma 
(QGP) that then during the expansion converts to a dense hadronic medium. 
Subsequently the chemical components (i.e. the abundances of the different 
hadrons) decouple and finally the system kinetically freezes out releasing the 
particles that reach the detectors. The early universe passed through various 
transition stages like the electro- weak transition, or the neutrino or photon 
decoupling, the latter leading to the micro-wave background radiation. 
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All such transitions share that they proceed during quite some time. The tran- 
sition can be accompanied by a strong change in the entropy density and thus 
releases quite some amount of latent heat, as e.g. during the QGP to hadron 
phase transition. Thereby the transition dynamics requires adjustments of the 
structural and thermodynamic properties of the system. 

Continuous transition processes can properly be described in the context of 
coupled transport equations, like the Boltzmann or e.g. chemical rate equa- 
tions, as e.g. presented in [2-7] . Still, sudden transition schemes are frequently 
employed, be it in the context of macro-dynamical descriptions or simply for 
a survey analyses of a wide body of the nuclear collision data. These concern 
coalescence pictures, which combine nucleons to composite nuclei at freeze- 
out [8] or coalesce quarks to hadrons in the deconfinement-confinement transi- 
tion, e.g. [9]. Such recipes do not only bypass the transition dynamics. In many 
cases they violate general principles as detailed balance, unitarity, conservation 
laws or entropy requirements [6,7]. Also the general freeze-out of the diluting 
system is mostly treated as a sudden transition happening at a suitably chosen 
three dimensional freeze-out hyper-surface in space-time [10]. Commonly one 
then defines a proper transition condition in terms of a transition temperature 
or transition density. In particular the latter freeze-out picture is widely used 
to analyse nuclear collision data in terms of thermal models. The achieved 
fits in temperature, chemical potentials and parametrised flow effects to the 
observed particle abundances and kinetic spectra are then frequently used as 
measured data that are taken as clues on the physics of the collision dynamics, 
cf. Refs. [11-14]. 

In this note we reinvestigate transition and decoupling processes under the 
perspective that these processes proceed during some considerable time span. 
Here we concentrate on the physics aspects of the continuous decoupling pro- 
cesses and shall defer the conceptual derivation, based on a flnal state interac- 
tion picture within the non-equilibrium formulation of quantum many-body 
theory as given by Kadanoff and Baym (KB) [15], to a forthcoming paper [16]. 

Starting from an exact microscopic formulation in terms of distorted waves 
[16-18] transparent semi-classical relations are obtained that describe a con- 
tinuous decoupling or freeze-out of strongly interaction particles. These for- 
mulae generalise the instantaneous freeze-out concepts, like the Cooper-Prye 
formulae [10], to finite decoupling time effects. Thereby the drain of energy- 
momentum and conserved currents from the source needs to be considered in 
order to come to an overall conserving scheme. The concept offers a direct gen- 
eralisation to two-particle coincidence rates known as the Hanburry-Brown- 
Twiss (HBT) effect [19] using distorted waves [20,21]. 

Earlier efforts [22-26] concentrated on introducing a finite depth of the decou- 
pling layer, across which a fiuid phase is coupled to a gas of frozen-out par- 
ticles. Those attempts simply parametrised the finite decoupling, rather than 
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inferring it microscopically [27] from the underlying decoupling processes, as 
will be considered here. The microscopic treatment permits to account for the 
individual properties of the particles, which e.g. leads to specific conclusions 
on the fate of short lived resonances during freeze-out. The latter is closely 
related to the generic balance between the local creation and the subsequent 
absorption during escape, which is finally responsible for an erosion of the 
memory on structural effects of the source. 

In the discussion section various decoupling phenomena, such as the chemical 
and thermal freeze-out, but also implications on phase transitions in nuclear 
collisions are reinspected under the particular perspective, that transition pro- 
cesses proceed during a considerable time span, across which the bulk prop- 
erties of the expanding system may significantly change. 

We use the convention h = c = 1 and formulate all relations for relativistic 
scalar bosons, for which standard covariant notation is used. 



2 Decoupling rates 



The invariant single-particle detector rate of a particle a (scalar boson) can 
be formulated as [16-18p] 



[27,) 
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Here the current-current correlation function n^'^'° encodes the property of the 
source, while the distorted single particle wave functions tpf^ describe the final 
state the propagation of the observed particle. This includes optical deflection 
and absorption effects through the real and imaginary parts of the retarded 
polarisation function n^(x,x') = iQ{t — t') ^ Jl{x), Ja{x') ^ . Expression 

(1) is rigorous, provided one knows the exact one-particle irreducible current- 
current correlation function and the corresponding exact final-state distorted 
wave ip^^l- The latter results from the outgoing scattering solution of the Klein- 
Gordon equation governed by Il^{x, x') with corresponding asymptotic bound- 
ary condition [28,29]. The technical and further conceptual developments will 
be the subject of a forthcoming paper [16]. 



^ The current field operator Ja{x) used in (1) formally results from the functional 
variation of the interaction Lagrangian with respect to the field operator 4>a{x) of 



particle a as Ja{x) = 5£™*(^)/5(/)^(x). 



3 



2.1 Semi- classical decoupling rates 



The semi-classical approximation provides a physically intuitive formulation in 
terms of classical paths. It either amounts to consider JWKB like approxima- 
tions for the distorted waves, leading to a covariant Hamilton- Jacobi problem 
for the classical action [30], or equivalently to perform the first order gradi- 
ent approximation to the Kadanoff-Baym (KB) equations [15,31,32]. In the 
small damping-width limit the approximate expressions for both cases are then 
unique and constructed through the bundle of covariant classical paths deter- 
mined by the real part of 11^. The paths obey the on-shell constraint, while 
the imaginary damping part of the classical action is treated perturbatively. 
Then each real path accounts for the entire local spectral width. 

In the following we will discuss two decoupling schemes: a) the local decoupling 
rate that can be even off-shell and b) the detector yield. For case a) the gain 
term, ns'^™(x,p)A(x,p), of the gradient expanded KB equation [15,31-34], 
where A is the spectral function, suggests the following local decoupling rate 



5^ = (2^n""°(^-^'^"(^-rtP=.-(-.rt- P) 

Here Pescape(a;,p) captures the probability that particles, created or scattered 
at space-time point x into a momentum p, can escape to infinity without 
further being absorbed by the loss part of the collision term. This formulation 
thus restricts the emission zone to the layer of the last interaction. Semi- 
classically in the small width limit V is given as 

nscape(a;,p)=e->^(^'P), where (4) 

oo 

X{x,p)= J V{x\p')dt\ and r(a;',p') = -ImH^(x',p')/Po> 

where the time integration is taken along the classical escape path starting at 
{x,p). The exponent x is also known as optical depth. At this point one may 
be tempted to assume that the local decoupling rate (3) together with (4) 
could also apply to the broad spectral-width case, upon generalising the clas- 
sical paths to appropriate off-mass shell paths. As yet a proof in compliance 
the overall conservation laws could so far not be giveif^ We'll nonetheless 
keep this intuitive formulation, which cares about both facets of damping, 
the attenuation of the flux as well as the corresponding spectral-width dy- 
namics, implicitly hoping that this form can be applied under much broader 
circumstances as used for its derivation. 

^ The occurrence of a back-flow term [31] objects to an exact formulation in terms 
of classical paths. A special approximation though was discussed in [31,33,34]. 
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Under the small width assumption that the spectral strength A(x,p) will be 
guided via the on-shell path, starting from {x,p) towards an on-shell detector 
momentum pa, leads to the following detector yield 



dKiPA) f d^xd^p ^^^^^ S^ipA-pAix^p^) 



d^PA J (27r 
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nrM„ Pescape ( ^ 1 5\p - p{x,Pa)). (5) 



Here Pa{x^p) denotes the corresponding mapping of the local momentum p to 
the detector momentum. The inverse mapping may neither be unique, 

nor even existing for some values of x (classical shadow regions which don't 
contribute) . The corresponding Jacobi determinant accounts for the focussing 
or defocusing of the classical paths due to deflections by the optical potential 
resulting from Re 11^. Owing to Liouville's theorem this Jacobian agrees with 
the van Vleck determinant [35] appearing as the pre-exponential factor in the 
JWKB ansatz for the distorted waves in (1). 

In the narrow width limit the above rates and yields naturally merges known 
continuous freeze-out pictures [22-26] of free gas kinetics, here however used 
microscopically with due account of the individual source terms and optical 
potential. The formulae encompass the completely opaque (i.e. strong inter- 
action) limit as well as that for weakly interacting i.e. penetrating probes. 
For the latter case final-state distortion effects are negligible and one exactly 
recovers the known Golden-Rule result 



dN ix v) 1 

~1^^¥^ = ^ nr"^(x, p)Al^ip) (for weakly interacting probes) (6) 
with Ar(p)=27r5(p2-m^), (7) 

here formulated in terms of the current-current correlator or polarisation func- 
tion, which still can account for non-perturbative multiple collision contribu- 
tions from the source [36]. 



3 Strong decoupling and freeze-out 



If one wants to describe the dynamics in terms of a two component scenario 
[22-26] with an interacting source and a decoupled (frozen-out) component, 
the above local decoupling rate (3) has to be supplemented by a corresponding 
description of the depopulation of the source. It causes a drain in particle 
number and energy and a recoil momentum. In fluid dynamical descriptions 
this leads to the following loss terms in the corresponding fluid cells 
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dNa{x,p) 



^^C.^-^Jl^pf'^. (8) 

These transfer rates result from the dissipative part of the gradient expanded 
KB equations, namely upon weighting the loss term with the charge Caa of 
particle a or with p'^, respectively. Here a labels a conserved charge. These 
recoupling terms have forms similar to those used in multi-fluid models, cf. [22- 
25,37-42]. Here, however, the decoupling of individual particles is addressed 
with their microscopic properties rather than an "anonymous" decoupling 
between two fluids. A complete treatment of the decoupling rate (3) together 
with the fluid drain terms (8) then provides an overall conserving scheme. Such 
recouplings lead to a gradual fading of the fluid phase, upon creating the freely 
streaming particle phase. Due to the a posteriori nature of the decoupling 
equations (3), since they involve the knowledge on the future, a solution of 
this coupled dynamics may only be obtained itcratively. As, however, pointed 
out in [27] such a two-component picture may not be quite consistent, since 
particles in both components are subjected to the same interaction dynamics 
and the tagging of the frozen-out particles is solely a matter of probabilities. 



Thermal source limit 



In thermal equilibrium the source function H^^™ can directly be expressed as 

Uf^(x,p) = -2 /th(x,p°) ImU^(x,p) = /th(x,p°) 2p° r„(x,p), (9) 

where ra{x,p) is the local damping width of particle a. This property leads to 
quite some compensating effect, which is frame independent. Namely, for large 
source extensions the integral over the F-dependent damping factors in (5), 
which define the visibility probability Pt, integrated along any path leading 
from the opaque interior to the outside, equates to unity 



r(t) e-^W = 1, where x{t) ^ j dt'T{t'). (10) 

This compensation is independent of the structural details of V and of the 
classical paths leading to the detector, along which the decoupled particles are 
accumulated. It is vahd as long as space-time changes in H or F are sufficiently 
smooth as to permit the underlying semi-classical considerations. In case of 
sharp changes in H the occurrence of reflection and diffraction effects [43] may 
partially spoil the picture. When rates drop smoothly in time, the visibility 
probability Ptit) achieves its maximum at 
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f(t) + r2(t) =0, where Pt(tmax) ~ r(t^,J/e. (11) 

Here the dot denotes the time derivative along the corresponding classical 
path. The corresponding decoupling duration At^ec approximately follows 
from the normalisation of the visibility function Pj to 

AW ^ j^^^ ^ ^yf-y (12) 

This defines a kind of decoupling uncertainty relation r(tmax) Atdec ~ e. 

Collecting all rates along the bundle of classical paths leading to the detector 
according to (5) provides 



(2^)^^^^!^ = / d'x 2/(i//th(x,p°) [§-,Aa{x,p) r,(x,p) e-^(-'^')] (13) 

_ . roc 

(14) 



where the square bracket expressions [. . .] are to be taken at the corresponding 
classical three-momentum p{x,pa) of the local detector orbit. This rate gener- 
ahses the Cooper-Frye formula [10] to a smooth decoupling situation dictated 
by the properties of the damping rate during the decoupling process. Similar 
continuous decoupling concepts, though based on the Boltzmann equation, 
were used in cosmology already for quite some time, see e.g. [44] and for a re- 
cent text book in this context [45], Chapt. 9. In (14) the space-time integration 
is expressed in terms of 3 dimensional hyper-surfaces layers a in the sense of 
curved coordinates. For each momentum each surface a can e.g. be defined as 
surfaces of constant r{x,p), while dx^ defines the world line direction towards 
the detector. 

In the opaque limit of the source and as long as the matter properties such as 
the temperature do not significantly change across the decoupling process the 
integral property (10) leads to an improved Cooper-Frye picture 



(2^) 



^dNaipA) 



d^Py 



CF 



d^a^p^dp' [^/th(x,pO)^f°(x,p) 



0"fo(RA) 



(15) 



with freeze-out hypersurface CTfoiPA)- As a special ingredient this result is 
compatible with Planck's radiation law for static (i.e. spatial) surfaces, which 
the Cooper-Frye formula is not. The tiny but essential difference is that the 
freeze-out surface does depend on the detector momentum pA [27,46] Thus 



^ Earlier attempts [24,47-49] tried to heuristically rescue this by a 0-function factor 
which is only non-zero for those particle momenta which permit the particles to leave 
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viewing e.g. the collision zone from different sides leads to different freeze- 
out zones. Therefore the mere assumption underlying the previous freeze-out 
pictures, that all particles touching a globally defined freeze-out hyper-surface 
do freeze out, is simply inappropriate especially for spatial surfaces. In order 
to recover a conserving scheme, though, the corresponding partial particle 
losses have to be accounted for by an appropriate surface rccoupling term 
for the source fluid, the latter to be derived from the drain rates (8). The 
Cooper-Frye-Planck formula (15) is written as to include the spectral function 
A^°{x,p) of the released particle right after freeze-out, which may be used 
energy differentially across the spectral width. For unstable particles, like e.g. 
vector mesons [50], it may then serve as an input for the subsequent decay 
e.g. competitively into dileptons and hadrons. 

The shining hyper-surface (Tfo(pA) in (15) can be determined for each observed 
momentum pA from the space-time points of highest brilliance as given by (11) 
to 



p^d,r{x{{t),p)+p' T\x{{t),p) 



= 0. 



(16) 



This provides a microscopic definition of a freeze-out criterion. A more general 
determination of the freeze-out surface, which refers to specific observables, is 
given in Eq. (28) in Sect. 5. 



4 Two particle final state correlations 



In straight generalisation of the microscopic definition of the freeze-out source 
(1) the two-particle source is given by the gain component of a double current- 
current correlation function of the two observed particles a and & [51] 



SabiXa, Va] Xb, Vb) = (4(^a )Mya)4i^b )Myb)) (17) 

with four-coordinates Xa etc. Accounting for the complete final state interac- 
tion (FSI) of both particles the coincidence yield expected at two detectors 
with momenta Pa and pb is then obtained by [17,20,21] 



hb{PA:PB) = I d'^Xad'^Vad'^Xbd'^yb'^yJp^iya, yb)Sab{Xa, Va, Xb, ?/b)^'t] (x^, Xb) 

^^ilf.)- (18) 

For details we refer to the review papers [51,52] and the historical perspectives 
[53,54]. 



the collision zone. This method lead to problems with respect to conservation laws. 
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In general the complete final state two-body waves ^ require a full scale three 
body problem to be solved. In this note we rather concentrate on the simplify- 
ing case, where the mutual interaction between the observed pair is negligible. 
Then the two terms in the ( ant i-) symmetric two particle wave functions still 
factorise and the result for the correlated yield over the single yields can be 
cast into a form 



C{pa,Pb) = 1 ± 





{a. 


■Qgain 




2 


{4a 


J^gain 




J^gain 





(19) 
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where now the single particle distorted waves ipp^ etc. enter. Here the matrix 
elements are most conveniently given in coordinate representation 



PA 



]-|-gain 



PB 



(iV(iVVii"n^)n^^"(^, r'Wplir'), etc 



PB^ 



(20) 



using the current-current correlator (1). In the plane wave limit of identical 
particles it simply recovers the standard result, cf. [51], valid for penetrating 
probes 



C{q,K) = l 



I o?V„ng^°(r„, Pa) I d^nUs^^'^in, Pb) 



PA=p]3=m? 



withg=pA-pB, K^-{pA+Ph)- 



(21) 
(22) 



For the strong coupling case again subtle compensation effect arising between 
source function and the attenuation along the escape path are expected for the 
matrix elements (20) in (19). In order to illustrate this we discuss the simplifi- 
cations emerging again in the semi-classical limit of large particle momenta and 
for small opening angles between the two detector momenta pA and ps- As Ja- 
cobian effects essentially cancel out in the ratio in (20) one can use the eikonal 
approximation evaluated along straight lines parallel to the averaged momen- 
tum K. Then the wave functions become '^^^{r) ~ e*^"''~2x(^'P») etc., where 
Pa and X the local four momentum and the corresponding optical depth. 
Using further the thermal limit property n^^'"(r,p) = 2p^fth{r,p^)r[r,p), the 
matrix element (20) becomes 



^ J d^fd'^K 2K°fth{r, K°) r(f, K) e^i^e'^^^^^) . (23) 

In the last step second order gradients were ignored in defining the averaged 
optical depth x(f, X). As a result the coincident yield can be described in some 



^PA 



■Qgain 
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quasi-free manner, i.e. using the plane waves form (21), thereby replacing the 
original source function n^^™ by an effective single particle source function 

S'eff(r,p) = ne-°(r,p)e->^('-'*'), (24) 

which accounts for the opaqueness resulting from the corresponding optical 
depths seen from the detectors. 

Alternatively to the analytic results given here one can sample the points of 
last interaction for the observed particles within transport simulations and this 
way define the proper emission zone profile in space and time, cf. e.g. [46,55]. 
In this case the simple yield ratio (21) is still applicable, modulo that one 
may have to account for the optical deflection of the particles due to the real 
part of the optical potentials, which maps the detector momenta to the local 
momenta. This opens the perspective for an imaging analysis of the active 
source region [56] or region of homogeneity [27, 46] . 



5 Analytic considerations 

The derived continuous decoupling relations raise quite a number of conceptual 
questions with respect to the observability of structural features of the evolving 
matter and the standard interpretation in terms of thermal model fits. First 
question to ask is, under which circumstances fingerprints of the structure of 
the matter are directly visible by certain probes? As all structural information 
is contained in the current-current correlation functions, there is a favourable 
answer: namely by penetrating probes. This has induced worldwide activities 
to observe photons and lepton pairs resulting from nuclear collisions, since 
these probes directly see the electromagnetic current-current correlator with 
its prominent feature to be influenced by vector mesons in the matter. 

For all strongly interacting probes, though, the interior is essentially opaque 
and one is left with surface effects, not only meant spatially but also in the 
context of time history. Then, however, the structural features are essentially 
washed out, since the structural parts of the decoupling rates encoded in 
the damping widths F essentially integrate to unity. This statement is also 
true in the near non-equilibrium case, since all non-equilibrium aspects can be 
encoded in a non-thermal distribution in (9), rather than /th- Thus fingerprints 
from incoherent decays of resonances with life times less than the decoupling 
time At dec may be less or differently visible. 

Dynamically expanding system have one robust feature, namely that on the 
mean all microscopic damping widths monotonically drop with density, ulti- 
mately proportional to the density for elastic scattering processes. This allows 
simple parametrisations for the time dependence of the microscopic damping 
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widths and this way leads to universal (scaling) behaviour between different 
physical systems. 

In the following we consider the emission of slow particles from a spherically 
expanding uniform fireball or from a Bjorken type expansion [57], once radial 
flow has built up. The emission is than essentially from a time-like hyper- 
surface. In order to come to simple analytic terms, which display the quali- 
tative features, we assume that the overall volume grows as V c<. in time, 
while the damping rates scale proportional to some power of the density p" 
times an averaged relative velocity. The case (n = 1) is essentially relevant 
for elastic scattering, while (n = 2) is more suited for inelastic processes. For 
simplicity the averaged relative velocity is taken proportional to the square 
root of the temperature as appropriate for massive particles, i.e. 

l/0C^^ {Vrel)ocT^/^ and T OC p"" {Vrel) oc t'" , (25) 

Tocl/-^''-^) and thus a = Sri + ^(k - 1). (26) 

Here we further employed some adiabatic relation between T and V with 
adiabatic index k = Cp/Cy- These assumptions yield damping rates at the 
decoupling peak and a corresponding spread in decoupling time, cf. Fig. 1, in 
the order of 

r(tdec) = «^ and Ataeo^- (27) 

For the examples presented in Fig. 1 we assumed some typical fireball values, 

decoupling probability [c/fm] 

0.15 
0.1 
0.05 

5 10 15 20 25 
t [fm/c] 

Fig. 1. Decoupling probability Pt{t) as a function of time for an expanding homo- 
geneous fireball for the three decoupling scenarios, where the rate depends linear 
or quadratically on the density (n = 1,2) or a mixed scenario as discussed in the 
text. The three curves slightly displaced to the right result from scaling all rates by 
a factor 1.5. 




11 



namely a freeze-out radius of -Rdcc ~ 6 fm and a collective velocity i?dec = 0.5 c 
for the nuclear collision case. Through (27) these collective parameters already 
determine the damping values at the decoupling peak, namely Tdcc = 0.5 c/fm 
100 MeV for inelastic processes (n = 2) and still 50 McV for scattering 
processes (n — 1). Beside these two options a mixed case with assumed 30% 
of n = 1 admixture at decoupling time tdcc = 12 fm/c might be appropriate 
for the kinetic or thermal decoupling, since in this case both, chemical and 
scattering processes, do contribute. Also here Fdcc emerges to about 100 MeV. 
Microscopically, of course, it is the time behaviour of r{t) that determines the 
decoupling conditions. 

The here obtained decoupling durations At^ec lie above 5 fm/c, during which, 
as a more universal result, the volume almost increases by an order of magni- 
tude. The robust microscopic figures in this respect, however, is the variation 
of the damping width between start (i) and stop (f) of the decoupling process 
(taken at full width half maximum). They vary by about a factor e^'^ ^ 1/15, 
i.e. Ti : Tdec : Tf fti es^ : 1 : e's^ ^ 390 MeV:100 MeV:26 MeV. The latter 
values of F representative for the above cases by themselves are tremendous 
news: during the entire decoupling process the resulting decoupling rates (or 
damping widths) are significant and of same order or even larger than the 
mean kinetic energies of the constituents. 

In order to see the sensitivity on the absolute damping strength a second set 
of curves is supplemented in Fig. 1, where all damping rates are scaled up by 
50%, correspondingly leading to slightly later decouplings. The corresponding 
curves scale to one another, if plotted in dimensionless variables and thus 
equally apply to other physical processes like e.g. the photon decoupling during 
the early universe. The latter is reported [45] to occur in a time window with 
the cosmological red shift factor varying between Z = 1300 and 800 with peak 
at Z ^ 1050, well in line with the time profiles deduced here. In the early 
universe case, however, the implied temperature drop during decoupling is of 
minutes harm for the nowadays observed cosmological microwave background 
radiation (CMB), since it is perfectly accompanied and exactly compensated 
by the cosmological red shift. Thus fluctuations in the CMB temperatures 
observed in different celestial directions result as tiny as on the 10~^ level. 

For the nuclear collision case there is a priory no such compensation to be 
expected. Rather from the above results one expects a significant spread in 
the thermodynamic conditions during the continuous decoupling or freeze-out. 
The time profiles of Fig. 1 can accordingly be transformed to the resulting 
probability distribution in temperature via 

Pt{T) = PMT))^, (28) 
where T{t) and its inverse t{T) describe the temperature change with time due 
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temperature distribution [1/MeV] temperature distribution [1/MeV] 



K = 1.5 



n 



= 2 



K = 4/3 



0.01 




50 100 150 200 250 300 50 100 150 200 250 300 



Fig. 2. Temperature distributions Pt{T) resulting from the three different decou- 
pling scenarios shown in Fig. 1. Left and right panels for EoS with adiabatic ratios 
«; = 1.5 and 4/3, respectively. 

to the underlying equation of states (EoS). For two simple example EoS with 
adiabatic ratios n = Cp/Cy = 1-5 and 4/3, respectively the corresponding 
distribution in T are shown in the two panels of Fig. 2. While k = 4/3 corre- 
sponds to an ideal gas of massless particles, e.g. photons in the early universe, 
the former interpolates between this massless case and the fully massive ideal 
gas case (M ^ T) for which k = 5/3. For orientation we used a tempera- 
ture value of Tdoc = 160 MeV at the decoupling peak around tdec = 12 fm/c. 
One also notices that due to this transformation one does not only obtain 
sizable spreads in the temperature distributions, but also that the peaks in 
T may significantly be shifted downwards compared to the decoupling value 
of Tdec = 160 MeV. Thus peak positions in T differ from those in t, illustrat- 
ing that depending on the observable the optimal decoupling condition may 
change because of the finite duration of the decouphng. 

Since the coupling widths F significantly change with volume and temperature, 
and hence drop in time, their absolute value at some given thermodynamic 
condition is less important than the relative speed by which these rates drop. 
This can be seen by comparing the results shown in both figures for the n = 
1 or also the "mixed" scenario, which have more a moderate drop of the 
coupling, relative to the faster drop for n = 2. The resulting distributions in 
temperature are vastly different. One obtains a much more narrow distribution 
peaking slightly below the anticipated decoupling T^^c of 160 MeV for the fast 
decoupling, while the mixed scenario relevant for kinetic decoupling peaks at 
much lower T with much broader width. 

Plots similar to the ones in Fig. 1 were obtained already in various transport 
calculations like in [2] for kaon production at SIS, or in [3,4] for the production 
of mesons and baryons at the CERN SPS, or in [5] for collisions at RHIC 
with comparable or even longer freeze-out durations. Still, those results did 
not quite trigger an alert towards a critical reinspection of the ubiquitous 
application of instantaneous freeze-out concepts. 



T [MeV] 



T [MeV] 
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6 Decoupling phenomena 



In the following we discuss several physical phenomena in the light of the long 
decoupling times. 

Chemical Freeze-out 

Over years nuclear collision data taken at the Brookhaven AGS, the CERN 
SPS accelerators and the RHIC collider were analysed in terms of thermal 
freeze-out models. In particular the relative abundances of the various particle 
species delivered the conditions on the chemical freeze-out in terms of temper- 
ature Tchcm and baryon chemical potential fiB- These fits provided data based 
freeze-out marks in a QCD phase diagram of T versus /is [11-13]. Thereby the 
T'chem did not only saturate with increasing collision energy at a value around 
165 MeV. For SPS and RHIC events good x-squared fits lead to relatively 
sharply determined freeze-out temperatures Tchcm for each collision energy, 
despite the here diagnosed freeze-out durations, during which the overall vol- 
ume increases by almost an order of magnitude. For ordinary EoS T steadily 
drops with increasing volume and one would expect a resulting spread in T 
similar to those given in Fig. 2. If the narrow temperature fits consolidate 
(see however [58]), could it be a sign, that the EoS is not such trivial and in- 
deed contains ingredients which provide a halt in the temperature drop during 
chemical decoupling? This would be an interesting chance, that results from 
the considerations given here. Phase transitions with considerable latent-heat 
release, such as the QCD to hadron gas phase transition, may be a possible 
explanation in this respect. 

Phase transition scenarios 

Phase transition involve a restructuring of the matter. The corresponding tran- 
sition rates and durations essentially obey similar uncertainty relations as Eq. 
(12) Ftrans Attrans ~ Alougsidc e.g. the dcconfinement-confinement transi- 
tion is accompanied by a strong reduction of the entropy density, implying 
a considerable amount of latent heat released during the transition. Thereby 
for the discussion given here it is less important, whether the transition is of 
first or second order or even cross over, a distinction of less importance for 
finite systems during a finite transition time. Important is that across a nar- 
row window in temperature a certain amount of latent heat is released. Such 
transitions then lead to a sensitive slowing down in the temperature drop dur- 
ing the transition process, as e.g. shown in a chemical rate equation model 
suggested for the description of the QGP to the hadron phase transition [6,7]. 
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Irrespective of the specific assumptions employed, on thermo- and transport- 
dynamically sound grounds this model provided qualitative insight into the 
transition dynamics. The dynamical investigations within a first order phase- 
transition scenario [59] show that such transitions take a considerable amount 
of time (about 6-10 fm/c depending on system size and beam energy). During 
this transition the temperature is essentially kept constant due to the release 
of latent heat, much like if coupled to a thermostat, while the volume increases 
by one order of magnitude though, cf . Fig. 1 in [6] . In such transitions the rates 
governing the abundances of the species in both phases are highly non-linear 
with density or temperature [6,60]. In order to provide the proper equilibrium 
phase they have to be governed by corresponding driving potential, such as 
the chemical potentials of the species resulting from the underlying EoS. Such 
scenarios will indeed provide much more narrow distributions in the resulting 
decoupling temperatures for the chemical abundances than in the simple EoS 
cases discussed here, cf. Fig. 2, since one expects the chemical decoupling to 
happen right after the confinement transition. 

Thermal Freeze-out 

For the kinetic or thermal freeze-out several competing effects enter, such as 
thermal motion versus flow, latent heat effects and the influence of the optical 
potential. 

Quite some compensation occurs between the thermal random motion and 
the collective flow that builds up during the expansion, cf. the results from an 
exactly solvable model [27]. An early freeze-out with small flow and a large 
temperature and a later freeze-out with larger flow and smaller temperature 
essentially lead to the same observed momentum spectra of the particles. Thus 
nearly invisible for a single component, flow effects can be isolated in compar- 
ing the spectra of species with different masses, cf. e.g. the analysis in [61]. 
Thereby the fits are not unique and therefore do not immediately contradict 
with a wider spread in kinetic decoupling temperatures. 

The special effect of the role of the escape probabihty (4) can be studied 
in comparing the freeze-out time structure of weakly and strongly coupling 
probes of similar or even identical mass. While the weakly interacting probes 
are essentially created during the entire collision process with an enhanced 
production during the initial phase for subthreshold processes, for the strongly 
interacting probes the freeze-out phase is shifted to the very late stage of 
the reaction. Transport investigations for the subthreshold production of K'^ 
(weak) to K~ (strong) in [2], or comparing the Q~ (weak) to proton (strong) 
production in [3] nicely confirm the here developed pictures with freeze-out 
durations that are well in agreement with the presented behaviour in the 
context of Fig. 1. That even absorption processes of rare probes can be counter- 
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balanced by rare multi-particle production processes was substantiated e.g. for 
anti-baryon production yields [62,63]. 

Finger prints of short lived resonances 

With the freeze-out duration a new time scale enters, which divides reso- 
nances into long and short lived ones. While long lived resonances survive the 
decoupling with their vacuum spectral function and subsequent vacuum decay 
patters, resonances with comparable or much shorter live times are affected 
by the decoupling dynamics. The question then arises, which possible signals 
do survive the decoupling process, e.g. for pions created through the decay 
of the A-resonance. Interesting in this context are the Coulomb corrected tt"*" 
and 7r~ spectra [64] and tt" data [65] taken at the GSI SIS. They show a two 
slope behaviour with a steeper rise towards small pion energies. Explanations, 
which attributed this low energy enhancement to the incoherent decay of the 
Delta resonance, arc disfavoTircd in view of relation (10), since the damping 
dependent factors integrate to unity across the escape path leaving the opaque 
zone. The physical origin of this is, that the visible layers are different in depth, 
such that the enhanced resonance production is essentially compensated by 
the corresponding reduced mean free path. The net result ultimately leads to 
Planck's law of black-body radiation in the opaque limit. Yet, such resonances 
do influence the corresponding optical potential of the emitted particle as en- 
coded in n^. In lowest order virial expansion the latter is determined by the 
corresponding scattering phase-shifts [66]. For the (tt. A?", A) system this con- 
tribution to the optical potential leads to some steepening of the pion spectra 
at lower cm. energies [67,68]. 

Composite-particle formation 

Composite-particle formation processes at freeze-out can be endo- or exother- 
mic. The formation of composites, such as deuterons or heavier nuclei for ex- 
ample, freezes degrees of freedom, thus releasing latent heat, that has then to 
be recoupled to the source according to Eqs. (8). Except for chemical freeze-out 
models this recouphng is mostly not accounted for e.g. in coalescence models. 
A microscopic description of composite light nuclei formation was formulated 
in [69] and further developed in [70] . Thereby composite particles are dynam- 
ically formed by multi-particle processes which respect the conservation laws. 

A particular example case of endo- and exothermic emission was studied in 
a conserving surface evaporation model applied to the decay of a thermal 
quark-gluon plasma glob [71]. While the emission of pions cooled the source, 
the emission of baryons simultaneously respecting energy and baryon number 
conservation indeed heated the source for the employed EoS. An interesting 
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feature in the context addressed in this note is, that this combined action lead 
to a stabihsed source temperature during emission, much hke for azeotropic 
distillation processes. 

Interferometry imaging and the HBT puzzle 

Prior to the RHIC experiments a wealth of model calculations appeared, cf. the 
reviews [53,54], that predicted large values for the ratio of two special HBT 
"radii" i?out/-Rside [72] with values beyond 4 or even much larger. However 
the finally measured data [73,74] arrived at values, which are similar to those 
already found at the CERN SPS, cf. [75-77], namely between 0.9 and 1.2. The 
origin of this puzzle [78] had many reasons. It partly resided in the use of simple 
freeze-out concepts, some times in conjunction with an inappropriate HBT 
recipe for the source's mean square lifetime, namely {{vxtY) ~ -Rout ~ -^side- 
As it ignores space-time correlations {v^ denotes the particles' velocity in out 
direction), this relation is valid only under severe restrictions, which prevent 
its usage for strongly interacting probes, see the discussion in [51], around Eq. 
(3.22). Besides space-momentum and space-time correlations already built up 
by the source encoded in H^^™, further correlation arise from the here discussed 
restricted optical view due to the final state damping as included in 5'eff, 
cf. (24). Together with a kind of collective Hubble expansion of the source 
[46,53,55,79] strong positive x-t correlations are built up, which significantly 
reduce i?out despite the long freeze-out times, since [80] 

^out = ((^ - V:ctf) - {{x - Vxt)f , while Rl^, = (y^) , (29) 

with spatial coordinates x and y pointing in out and side direction. The above 
averages are meant to be taken in the sense of Eq. (23). Thereby the opti- 
cal attenuation included in S^s leads to a possible breaking of some of the 
symmetries inherent in the original source distribution H^^'" and to corre- 
sponding shifts in the centre of gravity of -S'efif relative to that of H^^^". The 
here discussed considerations are nicely confirmed by recent hybrid transport 
simulations that appeared during the revision stage of this paper [46,55,79]. 
These model calculations did not only find the here advocated long freeze-out 
durations (above 10 fm/c). At the same time they provided a solution to the 
HBT puzzle, as the events emerged compatible with i?out/-Rside close to unity. 

In order to account for the distortion effects a technically involved distorted 
wave formalism was developed in Ref. [21] and applied to the analysis of pion 
coincidence data at RHIC. This approach however, did not quite respect the 
tight relation between the properties of the source H^^'° and the corresponding 
optical potential H^ determining the wave functions. Both ingredients were 
parametrised independently, although they are to be obtained from the same 
current-current correlation function as shown here. The absorption effects can 
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even be considered in a simpler way through the semi-classical effective source 
function S'efr, cf. (24), to be used in the plane wave formula (21). 

7 Concluding remEirks 

Prom general perspectives decoupling processes in dynamically expanding sys- 
tems show quite a universal behaviour. They proceed during a considerable 
time span, during which the damping rates drop by more than an order of mag- 
nitude implying alongside a sizable increase of the system's volume. The here 
formulated microscopic relations generalise instantaneous recipes and provide 
a conserving scheme for transition and decoupling processes in compliance 
with transport- and thermodynamic constraints such as detailed balance, uni- 
tarity and entropy requirements. Although the physical results are generally 
obtainable in much more detail from numerical transport simulations, the an- 
alytic formulae given here present a didactic view on such processes. This 
includes simple pocket formulae, which among others provide a kind of un- 
certainty relation between the damping rate of the decoupling particle at the 
decoupling peak and the decoupling duration. For the decoupling of strongly 
interacting probes thereby a subtle compensation effect emerges which leads 
to a generic behaviour. Like for Planck's radiation law it essentially wipes out 
the memory on the microscopic structural properties of the source. Despite the 
complications due to the finite decoupling time the presented scheme opens 
the perspective for improved future treatments of fluid dynamic calculations. 

Por applications to nuclear collisions the here diagnosed long decoupling and 
freeze-out times are both a challenge but also a chance: a chance to map 
out the thermodynamic properties of the expanding collision zone during the 
freeze-out of various probes. An observation of a quite narrow distribution in 
temperature, for example, could point towards effects that significantly slow 
down the temperature drop during expansion, and this way provide hints to- 
wards the underlying equation of state or possible phase transition effects. 
Therefore this note is meant as a stimulation to reconsider the analyses of 
nuclear collisions data in the light of the results and discussions given here. 
Promising steps towards this goal were recently achieved by hybrid model cal- 
culations [46, 55] , where the entire decoupling stage is treated within kinetic 
transport. 
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